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We consider the Higgs boson decay processes and its production and provide a parameterisation 
tailored for testing models of new physics. The choice of a particular parameterisation depends on a 
non-obvious balance of quantity and quality of the available experimental data, envisaged purpose for 
the parameterisation and degree of model independence. At present only simple parameterisations 
with a limited number of fit parameters can be performed, but this situation will improve with 
the forthcoming experimental LHC data. It is therefore important that different approaches are 
considered and that the most detailed information is made available to allow testing the different 
aspects of the Higgs boson physics and the possible hints beyond the Standard Model. 



Since the discovery of the Higgs boson last year by both ATLAS and CMS collaboration, it is clear that 
the study of its couplings is likely to give valuable information on the nature of the electroweak symmetry 
breaking and eventually to probe New Physics scenarios. Indeed, many of the currently studied BSM (Beyond 
the Standard Model) models can generically lead to deviations in the Higgs couplings, allowing for a constraint 
complementary to the direct searches for extra particles. However, to achieve such a goal, one must first make 
contact between the theory side (that is, the large amount of BSM models) and the experimental one (which 
are the results presented by the collaborations), and this implies the choice of a parameterisation. 

From the experimental point of view it makes sense to just parameterise Higgs physics in terms of observed 
quantities such as branching ratios and cross-sections. This is for example the case of the parameterisation 
proposed in Ref. [T], where the relevant cross-sections and partial decay widths are multiplied by a suitable 
factor. The advantage of such an approach is its simple link to the experimentally measured quantities. On 
the other hand, with such a choice, correlations among the different parameters are not explicit, in particular 
between tree level and loop induced observables. For example, a modification of the couplings to W bosons 
and top will also affect the loop-level couplings for the Higgs production via the gluon channel or the Higgs 
decay into two photons. Instead, we propose thus in [2] a parameterisation where the contribution of loops of 
New Physics to the H —> gg and H — > 77 modes is explicitly disentangled from the modification of tree level 
couplings. 

We must point that out that many studies have also been carried out relying on an effective field theory 
approach ([3, 3-6; ). While this approach has the important feature of allowing a full treatment of the radiative 
corrections (see [7] ) , it often relies on a large set of parameters, for which the current experimental accuracy is 
still lacking a bit behind. A similar study has also been recently carried out in [S]. 



In order to use the latest data from ATLAS and CMS collaboration ([9, 10J), one has to define the compatibility 
of a given model with data. The first approach is to use a % 2 test based on the signal strength /i which is simply 
the cross-section for pp — >• H — > XX normalised to the Standard Model expectation. Thus the \ 2 reads 



where \x is the best fit reported by the experiment, fi(n) the prediction of the model on the parameter point k, 
a the uncertainty, and i runs on all subchannels of each experiment. Then this x 2 (k) is compared to a standard 
X 2 distribution with n degrees of freedom, where n is the number of subchannels, in order to determine if the 
model is excluded or not. However this method suffers from a few shortcomings : first the fii(n) does not 
correspond to the inclusive cross-section, but to the exclusive one. For instance the H — > 77 decay mode in 
ATLAS is divided into 11 subchannels, which amount to as many exclusive cross-sections. To compute exclusive 
cross-sections, one needs the experimental fractions per mode £?, or equivalently the efficiencies : 



I. HIGGS COUPLING PARAMETERISATIONS 



II. USING EXPERIMENTAL HIGGS RESULTS 
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where Cf represent the SM fraction of the production mode p among all production modes in the subchannel 
i, and [i xx the signal strength of that specific production mode in the final state XX. Although the (f are 
computed by the collaborations themselves, there are not always publicly available. Another inconvenient raised 
by eq. [T] is that in adding all subchannels together we implicitly assume that they are not correlated. While 
this is the case of statistical uncertainty it is certainly not true for systematical uncertainty, let alone theory 
uncertainty. 

A. Improved \ 2 method 

There exists nevertheless a way to improve the statistical test : indeed the collaborations have released 2D 
plots of the \ 2 in each decay mode, showing the % 2 of this decay mode as a function of Ut gg H/ttH! Pvbf/vh) ■ 
Here £t gg H/ttH stands for a common signal strength for both gg — > H and pp — > ttH and fJ-VBF/VH f° r both 
VBF and VH. By approximating the \ 2 to a gaussian, we can trade the simple \ 2 test to a new one, defined 
as 

/ \ / u xx - u xx \ 

-v 2/ (k\ - V I r, xx - i, xx r, xx - n xx \ v- 1 i ss H / ttH 'ggHAtH \ (o\ 

X \ K J — / , \ MggH/ttH ^ggH/ttH' ^VBF/VH FVBF/VHJ v \~XX XX j ' 

XX WVBF/VH ~ ^VBF/VH' 

where V~ x is the inverse of the covariance matrix, deduced from the experimental plots. The advantages of 
the new method are straightforward : one does not need any knowledge of the fractions £f and moreover, most 
of the correlations between production modes are taken into account. 



B. Shortcomings of the improved method 

It is however clear that in order to use the improved \ 2 one has to make some assumptions. Indeed we are 
collecting 4 productions modes (gg — » H, VBF, VH and ttH) into 2 parameters (/% g H/ttH,Mv bf /v h) ■ However 
the requirement is not so stringent since any model abiding by custodial symmetry will feature an identical 
rescaling for both VBF and VH cross-sections. Concerning gg — > H and ttH , so far most channels are not 
sensitive to the latter because of its small cross-section, so assuming an identical rescaling does not affect much 
the prediction. 

Another issue is that the \ 2 reported by experiment is not a gaussian, and though this approximation may 
be accurate near the best fit, it goes worse as one moves away from it. However since there is no analytic form 
of the real \ 2 ''■> there is no easy alternative to the gaussian approximation. 

C. Experimental input 

The experimental input consists in the best fits and covariance matrix for 77, WW and tt decay modes 
for CMS and best fits and uncertainties for all subchannels in the remaining decay modes of CMS and all of 
ATLAS. The reason why we did not use improved % 2 for each decay mode is first that for decay mode ZZ and 
bb there is no much gain since the former is an exclusive channel and the second relies entirely on associated 
vector boson production. Second, all data from ATLAS was not available at the time to carry out the full 
improved \ 2 method. Those data were extracted [51 110] . and results for improved x 2 are shown in Table [l] 
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FIG. 1: CMS results in H — > 77, H — > WW and H — > tt channels. 
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III. 



PARAMETERISATION 



Our first parameterisation will be tailored to BSM models which alter mostly the Higgs physics via loop 
effects. This is a generic feature of models where there is no much mixing between the SM Higgs and any other 
scalars, but which feature extra particles light enough (light superpartners in Supersymmetry, light vector-like 
fermions in extra dimensional theories, and so on). In which case the effect of in each BSM set-up can be 
parameterised by the pair (k 9 , k 7 ), that we define as the amplitude of new particles contributing to the partial 
widths T gg and T 77 normalised by the SM top amplitude : 



. 77 



99 
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A t (T t )(l + Kj) + 
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(4) 
(5) 



where dots stands for the contribution of light quarks, Ax(tx) are usual SM amplitudes and contains the 
QCD NLO corrections (see [5] for details). 



A. Results 



For reference, sample points for the following models are indicated: 

- [♦] fourth generation where the result is independent on the masses and Yukawa couplings. 

- [*] Littlest Higgs [IT], where the result scales with the symmetry breaking scale /, set here to / = 500 
GeV for a model with T-parity. 

- [a] Simplest Little Higgs [12] . where the result scales with the W mass, also set here to rayfi = 500 GeV 
for a model with T-parity; 

- [■] colour octet model [13], where the result is inversely proportional to the mass mj = 750 GeV in the 

example; 

- [(g)] 5D Universal Extra Dimension model [2], where only the top and W resonances contribute and the 
result scales with the size of the extra dimension (here we set tukk = 500 GeV); 

- [^k] 6D UED model on the Real Projective Plane [15], with tukk = 600 GeV is set to the LHC bound [16] : 

- [•] the Minimal Composite Higgs [T7] (Gauge Higgs unification in warped space) with the IR brane at 
1/R' = 1 TeV, where only W and top towers contribute significantly and the point only depends on the 
overall scale of the KK masses; 

- [▼] a flat (W at 2 TeV) and [6] warped (1/R' at 1 TeV) version of brane Higgs models. The result only 
depends on the overall scale of the KK masses. 

One must note that each model will not be represented as a point in the (k 9 ,k 7 ) plane, but rather by a 
line starting at the SM point (0,0), since they all have a decoupling limit, except for the 4 th generation. We 
show in figure [2] the one and two sigma excluded regions, and the position of the models with respect to those 
exclusions. As one can see, the 4 th generation lies well away from the compatible region, and so do some of the 
benchmark of the other models. In particular, the 6D UED benchmark that we used was chosen so that the 
heavy scale was at the limit of the direct searches for extra particles, and we see that in this case, the indirect 
bounds form Higgs physics does much better than the direct search. 



IV. CONCLUSION 



We have shown how to go beyond the simplest methods when using experimental data to constrain Higgs 
couplings and account for part of the correlations between measurements. We have presented a parameterisation 
([HHH]) and we showed how it can be used for testing and putting exclusion limits on models of new physics 



Toyama International Workshop on Higgs as a Probe of New Physics 2013, 13-16, February, 2013 



4 




FIG. 2: «;' and K gg at the LHC for a Higgs boson with mjj = 125 GeV. The two solid lines correspond to the SM values of the 
inclusive 77 channel (A), and the vector boson fusion production channel (B). On the left panel, the fit using ATLAS data. On 
the right, the fit using CMS data. Both fits use 77, ZZ and bb channels. Darker (lighter) blue are the 1, 2 a limits. 



beyond the Standard Model. In particular our parameterisation is tailored to investigate BSM models, keeping 
track of the specific correlations among the parameters. It also allows more easily to interpret mass limits and 
contributions to the loops giving the effective Higgs boson vertices. 

We also performed 2 parameter fits of the CMS and ATLAS results using all available channels, showing that 
they already include all the necessary information and are therefore a good approximation at this stage. More 
precise measurements of extra channels will require the inclusion of more effective parameters. 
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